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A PERFECT -GAS ANALYSIS OF THE EXPANSION TUNNEL, 

A MODIFICATION TO THE EXPANSION TUBE 

By Robert L. Trimpi and Linwood B. Callis 
Langley Research Center 


SUMMARY 


A perfect-gas analysis is presented for an apparatus consisting of a basic 
expansion tube at the downstream end of which a nozzle has been added. The 
resultant apparatus, named the expansion tunnel, is shown to have the following 
advantages when compared with the basic expansion tube: increased testing 

time, larger initial test-gas slug length, higher efficiency, and reduced 
secondary-diaphragm bursting problems. Principal disadvantages are require- 
ment of an additional nozzle, and the requirement of a larger ratio between 
maximum and minimum pressures in an operating cycle. Authors conclude inherent 
advantages more than compensate for disadvantages. 


INTRODUCTION 


For the past few years scientists at the Langley Research Center have been 
investigating various modifications to the basic expansion tube described in 
reference 1. Effort has been concentrated on those variations which were more 
directly aimed at alleviation of the anticipated principal undesirable features 
of the expansion tube; namely, the short initial length of the test-gas slug 
before diaphragm rupture, the bursting of the secondary diaphragm, and the 
short test time. Both experimental and theoretical investigations have been 
conducted, and the latter include considerations of several modifications for 
both real and perfect gases. Most of the experimental results obtained to date 
in the pilot expansion tube at the Langley Research Center (unpublished) are 
for expansion- tube operation; these results appear encouraging. 

This report will describe the perfect-gas analysis of the configuration 
that the authors believe holds the most promise for reducing the aforementioned 
drawbacks. This configuration, mentioned in references 1 and 2, is called an 
expansion tunnel and consists of a basic expansion tube to which at the down- 
stream end a nozzle has been added. (See fig. 1.) Thus the test fluid is 
processed first by an unsteady expansion in the accelerating chamber and then 
by a steady expansion in the nozzle. A group at the Von Karman Laboratory of 
the Arnold Engineering Development Center has also been investigating both 
theoretically and experimentally a different configuration which has a nozzle 
located just after the secondary diaphragm. For such a configuration the fluid 


is processed first by a steady expansion and subsequently by an unsteady expan- 
sion if the accelerating chamber is several orders of magnitude longer than the 
nozzle. If these length restrictions are not satisfied, the test gas is proc- 
essed simultaneously by both steady and unsteady expansion waves with the 
result that the gas state at the test section continually varies with time. 

The relative prior occurrence of the unsteady expansion as compared with the 
steady expansion has an extremely important bearing on the subsequent charac- 
teristics of the apparatus. Consequently, even though both modifications 
contain nozzles, their operation and performances are very different. 

The analysis herein is restricted to the perfect-gas assumption which per- 
mits the pertinent equations to be expressed in closed form. Such equations 
are valuable since important trends and influences can often be simply 
extracted and examined critically. A somewhat parallel real-gas analysis 
(unpublished to date) has also been executed. The real-gas analysis verifies 
the perfect-gas trends of this paper although the magnitudes of the variations 
naturally are not identical. 

For the convenience of the reader an index to the figures is presented as 
table I . 
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quasi-one -dimensional-flow area 


geometric nozzle area ratio measured normal (perpendicular) to 
axis 

speed of sound 


for linear area nozzle and conical nozzle, respectively 


driver-gas specific heat at constant volume 
test-section (nozzle exit) diameter 
accelerating-chamber diameter 
nozzle entrance diameter 
driver-chamber energy (eq. (72)) 
driver-chamber energy parameter 
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gas state before nozzle (fig. l) 

ratio of nozzle losses for linear radius and linear area 
nozzles (eq. (62)) 

gas state in test section 

parameters defined in eqs. (79) 

gas state (fig. 10) 
total enthalpy 

functions defined in eqs. (43) and (56) 


initial state of driver gas before arc discharge 
constants defined in eqs. (44) and (57) 


Z N^S r 


lengths of driver, intermediate chamber, accelerating 
chamber, nozzle, and dump-tank sections of expansion 
tunnel (see fig. l) 


length of expansion chamber in shock tunnel 


flow Mach number 
test-section Mach number 


primary-shock Mach number, 

points on wave diagram (fig. l) 
integer 

Riemann parameter (eq. (9)) 
static pressure 

Reynolds number based on test- section and accelerating-chamber 
diameters 

absolute temperature 
time 

time increments (see figs. 1 and 13) 
velocity of primary shock wave 
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U Sll 

velocity of secondary shock wave 

c? 

H 

O 

velocity of tertiary shock wave 

Uo 

b R 

velocity of reflected shock wave 

u 

flow velocity 

X 

distance in flow direction 

z 

initial length of test gas slug in intermediate chamber 
(eq.. (67)) 
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fi(Mf) 


r^M) =1 + 21 

S 
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2 

7 

ratio of specific heats 

S 

boundary-layer thickness 

■n 

overall efficiency > r )o T U 1 

IT a 


ideal efficiency 


nozzle time loss efficiency 

^d ' 

nozzle capture efficiency 

e 

flow angle 

e w 

wall angle 

M- 

viscosity 

1 

nozzle coordinate in flow direction 

P 

density 
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ft(M) = r x (M) [r 2 (M)] 


- 1/2 


CD 


exponent in viscosity-temperature relation 
Additional remarks regarding notation: 


© 

( )i 

( ) 


ETun- 


( ) 


NRS 

( )* 


denotes gas in state i in cycle 

subscript signifies quantity is to be evaluated in ith state 

( ) ET refers to expansion tunnel , expansion tube, or nonreflected 

shock tunnel, respectively 

refers to critical conditions when At^- es -(- = 0 (that is, 
nozzle time losses equal to Atg) 


THEORY 


A description of the expansion- tunnel components and operation is more 
easily explained with reference to figure 1 which shows a schematic of the 
tunnel and a distance-time or wave diagram of the operating cycle. The tunnel 
has five basic sections: driver chamber, intermediate chamber, accelerating 

chamber, nozzle, and dump tank. Three diaphragms separate the first four pre- 
ceding sections. Thus, the expansion tunnel is an expansion tube to which has 
been added another diaphragm and a nozzle. The test fluid, which is initially 
in the intermediate chamber, is processed first from state © to state (2) by 
the primary shock wave, next to the state (e) by the unsteady upstream expan- 
sion wave, and finally to the test state (f) by passage through the nozzle. 

A perfect -gas theory is developed to illustrate the important facets of 
the expansion tunnel. The assumption of strong shock waves is used for all 
shock waves. Approximate equations are also often given in terms of test- 
section Mach number and area ratio for the asymptotic limiting cases of large 
nozzle exit and entrance Mach numbers. Each phase of the cycle is treated by 
considering in turn the processes of the gas initially in the intermediate 
chamber, accelerating chamber, driver chamber, and so forth. Readers not 
interested in details of the derivations, flow processes, and so forth, can 
proceed to the section on "Discussion" wherein the main points are considered. 


Test-Gas Processes 


The test-gas processes are found by working backwards from the desired 
test- section conditions © through the intermediate states © and (2) to the 
initial charging state (T) . The test-gas value of y is implied when y is 
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used without a subscript. All curves shown in this report are for a test gas 
with y = 1.4. 

The states © and (?) are related by the familiar isentropic quasi-one- 
dimensional steady-flow equations (see ref. 3)' 
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(3) 

( 4 ) 


The variation with test- section Mach number M-^ of the nondimen sional 

parameters on the left-hand side of equations (l) to (4) for the perfect test 
gas (7 = 1.4) is given in figures 2(a) to 2(d). The area ratio A varies from 
the basic expansion tube value of unity to 1000. The curves are terminated by 
a short dashed line at the lower Mach number end when the expansion fan has com- 
pletely vanished (M e = Mg = 1.89) and the apparatus is then operating as a non- 

reflected shock tunnel. (Of course, in a nonreflected tunnel one would test in 
gas initially in the accelerating chamber rather than in the intermediate 
chamber. ) 


The asymptotic approximations for equations (l) to (4) are obtained by 


assuming 



and 




Then from equation (l) 
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and from equations (2), (3), 


and (4) 
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Pf 


These approximations are indicated at the higher values of Mp in fig- 


ure 2. Note (fig. 2(h)) that for Mp ^ 20, the exact values of — ^ 0-9; 

u f 

therefore, only a small percentage increase in velocity occurs in the nozzle. 


Since the process from © -» © -» © is simply the expansion tube 
process of reference 1, the equations of that reference are applicable to the 
test gas by substitution of state 0 for state (^) • However, a brief outline 


of the necessary steps is as follows: 


The strong shock assumptions 


-£■ » 1 


require that 
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(6) 


m 2 


2 

7(7 - 1) 


(7) 


*%q 



7+1 ^2 ^ 7+1 ^2 
2 7 Pq 2 aq 


(8) 


(Eqs. (5) to (8) are identical to eqs. ( 3 ) to (6) of ref. 1.) The value of 
M 2 is approximately I .89 for 7 = 1.4. The Riemann parameter P is con- 
stant across the upstream wave (eq. ( 7 ) of ref. l) ; thus, 

P = — a 2 + u 2 = a e + u e (9) 
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Consequently, 



P2 
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Combining equations (2) and (4) with equations (lO) and (ll) yields 


where 



a 2 a 2 a e 
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Other pertinent parameters in state (§) are then 


Ug Mg ag Mg 

~ Mf *f ~ % rq(Mg) 


% = He = 

n(Mg) 

l 2 

n(M 2 ) 

Hg Hg 

fi(Me) 


fi(Mf) 


( 12 ) 


(13) 


(14) 


(15) 


(16) 


Plots of equations (12), (l3)> (15 ) > and (l6) appear in figure 3- The 
values of ag, Ug, Pg, and Hg all are increased as A is increased for a 
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fixed 


Mf . An approximate expression for 3 with the restrictions 


Mf 2 » 




M e 2 » 


- 1' 


and for Mg of the order of 


7-1 


is 





X M f)L 


7-1 


Mf 



0 


1 + 


7-1 

(A) 2 - 1 

r l( M f) 


(17a) 


(17b) 


(17c) 


Both exact equations (l4) and approximate equation (17) values for 0 are 
plotted in figures 4(a) and 4(h) . The value of 0 is increased by increasing 

A or decreasing Mf and approaches unity only for small A. However, at 
M = 50 and A = 1C)3, (3 is reduced to approximately 1^. The importance of 


the parameter 0 lies in the fact that it directly relates the requirements of 
the expansion tunnel to the expansion tube for a given Mf. (Note that in this 

report Mg ~ I.89 is constant for all apparatus.) 


The amount that 0 exceeds unity is an indication of the loss when the 
final expansion process is steady (expansion tunnel) rather than completely 
unsteady (expansion tube) . As a consequence of the reduced-enthalpy multipli- 

/% l\ 

cation in the unsteady expansion — « — |, it will be shown that not only must 

\h 2 

_ 2 _ 

the initial charging pressure be increased pp °c 0 but also the primary- 
shock Mach number ^Mgp a p^ . 

The primary-shock Mach number can be evaluated by combining equations (8), 
(12), and (15); thus. 
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r l(Mf) 
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(18) 
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The initial charging pressure ratio P^/Pf I s easily obtained by appli- 
cation of the gas law and the limiting strong shock density ratio 


2 


Pi 

Pf 


/a f \ 2 _ Pq /Pq\/P2\ _ 7 - 1 

1 ’ I ( Mf ) « 

\ 1 J Pf ylpf ) 7 1 

C 1 

^1 


7-1 


(19) 


Figures 5 and 6 are plots of equations (l8) and (19). Figure 5 indicates 
almost a constant increase in with A independent of Mp, whereas 

figure 6 shows a marked dependence of Pq^Pp on Mp, the relative penalty 

associated with increasing A declining with larger Mp. 


Accelerating-Gas Processes 

The cases of accelerating gases with specific heat ratios of both 

1.4 and 1.67 are considered. The initial charging pressure p-^ is first 
determined by application of the strong shock approximation for P 21 /P 11 
together with equation (4) and the equalities p g = p 2 q and u e = u 2 p: 


p 21 

P 11 


7 lll 

(Ui - 1 ) 

/M 

2 

\ u f j 


Mp C 



( 20 ) 


27 

P ll/ f \ _ 2 ( u f\ 1 / a e\ 

p f \ a ll/ 7 ll( 7 ll + 1 ) \ u e / Mp 2 ^!/ 


( 21 ) 


An 

re suit s 


evaluation of the absolute pressure level as influenced by and 

from manipulation of equation (2l) into the form of equation (21a) : 
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27 
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(21a) 


and in approximate form 
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(21b) 


P ll/ a f \ _ _2_ 
P f \ a i 


7i 7 1X + 1 


1 pj /Molecular weight of gas in state 

molecular weight of gas in state (O^ 


. ]_44 

In figure 7 > equation (21a) is plotted for values of - — = and 

7 1 n = 1.67- This value of a qq/ a q is representative for T 1 = *11 with per- 
fect air in region © and helium in region (l^ . Since for many testing 

Pu / a f \ 

purposes ap and might be nearly equal, the ratio I xs a 

direct measure of the charging pressure. The pressure level in region 
increases rapidly with A (eq. 2l(b) and fig. 7) and reduces the low vacuum 
pumping requirements for the intermediate chamber. If air were also used in 

p ll/ a f ) 2 


region (11) , the values of — 1 would then be 0.154 times those of 

W p f \ a l/ 

figure 7- 

The secondary-shock speed follows from equation (8) : 
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Uf f 


( 22 ) 


U c 


7u + 1 


>11 


u., 


(22a) 


Dump-Tank Processes 

The conditions in the dump tank are found by assuming a perfect nozzle 
start and that the dump tank has an area equal to that of the nozzle exit. 

Such a start requires that there be no upstream compressions generated by the 
passage of the secondary shock Ug through the nozzle. (See ref. 4.) Thus, 

the state (lo) must be selected so that after acceleration to the state (20) by 


the tertiary shock Ug , the following conditions are satisfied: Pg~ = p 


3 10 

u 20 = u 30‘ State (30) is defined by the steady-flow nozzle expansion of 


J 30 J 


state (23) . The relations between states (21) and (30) are analogous to those 

/ 

between states (e) and © except that Mgp is constant at 


'Ill'll " X ) 


11 



Thus with the strong shock assumption for 
nozzle charge pressure is 


Ug , the expression for the initial 


PlO = /£loV£3oVPe\ 
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The approximate values of 


u 30/ u 21 and p 30/ p 21 


are 
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(26) 


Substitution of equations (3a ) , (4a), and (24) to (26) into equations (23) 
results in 
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P T 0 /a f \ 2 

The nozzle charging pressure parameter — — [ ] of equations (23) is 

p f \ a 10 / 

plotted in figure 8 for 7 =1.4 with 7 -^q and equal to 1.4 and I. 67 . 

As indicated by equation (23a), this pressure parameter is only a weak function 
of area ratio when 7 = 7 and A » 1 . (See figs. 8(a) and 8(b).) Lower 

charging pressures are required for larger values of A when 7^7 


11 


A parameter illustrative of the absolute pressure level for the case of 
air as the test gas with helium or air for the other charging gases would be 


p io/ a f 
p f \ a l 


, which is found as 


p io/ a f \ 2 = p io/ a f \ 2 / a io N 


p f \ a 


\ a iOy 


( 27 ) 


This parameter is plotted in figure 9 f° r the ratios 


*10 

*1 


=1.0 for 


a i n i44 , , 

7 n „ and = for 7 , _ = 1.667- These sonic-speed ratios are appro- 

xu a-^ 49 xu 


priate for T 1 = T 10 with air or helium as the nozzle gas in state (10) . The 

lower pressures (higher vacuums) are required for air in both states (10) 

and . The combination of helium in state ^0) and air in state per- 

mits higher initial p 1Q (for Mf > 25, A > 25) than if helium were used in 

both states (10) and (G) . The case of 7 -^q = ^ and 7 -^ = ^ requires a 

p-^Q reduction of about one order of magnitude between A = 1 and A = 100, 


and the value of p. 


10 


5 X 10 -if p£ at 


10 . 


The exact expression for the tertiary shock velocity U q (eq. ( 28 )) may 

D 10 

be put into a simpler approximate form by using the strong shock relation 
together with the assumption of a perfect initial nozzle start and equa- 
tions (3a) and (24): 
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When the test gas flows through the nozzle from state © to state .©, 
the velocity u f ^ u^q and p f / in fact, from equations (3a), ( 4 a), 

( 2 k), and (25) 
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Since for A 1 , 


Hjo 5 M 2 ! 


'll (’ll 
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it follows from equa- 


tion (29) that u^q > Uf. Also, for 7 H 7^, p^ 0 < p f from equation ( 30 ). 

Consequently, the unsteady wave system generated by the passage of the test- 
gas — accelerating-gas interface ^between state © and state (2^) through the 
nozzle must attempt to adjust the states © and (30) to a common velocity and 
pressure. A schematic sketch and distance- time (wave) diagram of the processes 
leading to the establishment of the states (f) and (30) are shown in fig- 
ures 10 (a) and 10 (b) . For a perfect initial nozzle start, state (3^) is iden- 
tical to state ( 4 ^) . The pressure and velocity of states (f) and (§d) 

(fig. 10(c)) must now be matched at a common point (g,50) by means of an 


upstream wave, state (f) -e> state 


state 


state 


with > 0. 
du 


with — < 0, and a downstream wave, 
du 

These conditions are indicated for 


cases I, II, and III by the dashed lines in figure 10 (c). If both adjusting 
waves are isentropic, the equation governing the wave strength Pg/ p f 


is 
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Equation (28) is only dependent on 7, Mp, and A. Solutions of 

equation (3l) showed that, in general, cases II and III of figure 10(c) applied 

5 

with all reflections as expansion waves. Only for the condition y 11 = — and 

3 

M 0 -> Mg did the wave system of case I occur. Since for this condition the 
maximum pressure ratio P^q/^O was onP y 1-6 f° r A = 10^ and Mg = Mg, the 
isentropic equation (3l) is still a close approximation even for case I. 

The wave strengths Pg/Pp ^ or 7 = 7]_p = 1-4 are plotted in figure ll(a) 

and for 7 = 1.4, 7,, = — in figure ll(h) . Note that ^ ^ 1 so that 

xx J Pp 

upstream expansion waves are generated which will not disturb the expansion 
flow at the nozzle-exit model location. There are no solutions in figure ll(b) 
for the larger area ratios because for these cases no positive values of pp 

will satisfy equation (28); that is, case III of figure 10(c) falls below the 
p-axis. Physical considerations naturally limit p g and p™ to zero; there- 
fore, the wave system for these cases creates a perfect vacuum which grows in 
extent with time since u g < u=>q at p = 0. This high- vacuum-producing mecha- 
nism may offer interesting possibilities for other experiments. 


One more wave interaction must be considered in the dump tank, namely the 

end reflection of the shock wave with a velocity of U Q . If it is assumed 

b 10 _ 


that the wave is unaffected by any overtaking waves and that the state 


is 


the state into which the reflected shock wave with velocity Ug-^ will advance, 
the following equation applies: 
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If the strong shock approximations are used, equation ( 32 ) simplifies to 
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Also under the strong- shock assumption, the pressure after the reflected 
shock p^g is 


P RS = £ 7 10 ~ 1 
P 20 7 10 “ 1 


(33) 


P RS = %3 £20 £l0 
P f ’ P 20 P 10 Pp 


(34) 


Substituting approximate equations (23a) and ( 28 a) into equation (34) 
yields the large-area-ratio asymptotic expression 
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Consequently, p R ^ is of the same order as p.p and will be (for large 
area ratios) essentially independent of A if 7 -^ = 7 , whereas p R g will 
decrease with increasing A for 7 ,-^ > 7 . Solutions of the exact equation (34) 
are not shown but follow the trend described for large Mp. However, as Mp 
decreases, Ppg/Pf increases somewhat. 

Driver-Gas Processes 


p 4/ 


The driver-gas pressure level is evaluated by first finding the ratio 


p^ across the driver-gas unsteady flow expansion: 
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The ratio of driver-gas pressure to test-gas pressure results from combining 
equations (12) , (13), and (35): 
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The ratio of driver-gas pressures required for expansion-tunnel operation 
to that for expansion-tube operation = lj , for the case of identical test- 
section conditions, is then 
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27 

fi - • 

1 H 

jro. 

1 

7-1 

2 

\ a 4 /ET 


1 74 " 



2 

' wETun 


(37) 


(37a) 


u p 

If a driver gas with a very high speed of sound is used so that — « 1, 

a 4 

the expression in brackets approaches unity. If, in addition, the approxima- 
tion (l7c) is used for 0, the driver-gas pressure ratio becomes 


( P4 )ETun _ n7 -l 

(Net “ 


( P4 )ETun 

( p ^)et 


7-1 

1 + _£Z — ( A ) ..., 2 + 

7 - 1 T^Mf) 
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+ 

(A) 2 - 1 

1 

H 

ru 

r l( M f) 


(38) 


(38a) 
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If the assumption is added that Mf » 


7 - 1' 


N 


7-1 


7-1 

•rx 2 


— 1 2 


N: 


ETun w x + 4 7 (A) - 1 + 4y(y + l) [(A) - 1 

ET (7 - l ) 2 ^ (7 - l)^ fif . 2 


+ • 


(38b) 


Driver-gas pressure ratios from equation ( 36 ) are plotted in figure 12 for 


5 a 4 144 / t 4 t 4 

the case of 7. = =-, — = / — , with values of — = 10. 25, and 50. The 

4 y a f 49 f T f T f 


increase in driver-gas pressure with increasing nozzle area ratios is evident 
together with a decrease in p^ as a^T^) increases. Such a trend was 

expected from shock-tunnel experience. The decreasing difference between 
tunnel and tube pressures as Mf increases is also apparent, a trend predicted 

simply by the approximation equations ( 38 ). At large Mf, ^ P 4 /Pf ^ l° r A = 10^ 

is about ten times ^P4/Pf^ for A = 1. 


LENGTHS OF COMPONENT SECTIONS AND TESTING TIMES 


Accelerating Chamber 

The accelerating chamber length Iq 2 is the fundamental one in the 
apparatus since it determines the testing time. Let At 2 be defined as the 

time interval at the nozzle entrance between the arrival of the test gas and 
the trailing edge of the upstream expansion fan; that is, it is the time for 
which the state © exists at the nozzle entrance (figs. 1 and 13 ) and would be 
the test time for a zero-length nozzle. The nondimensional equation for At 2 
is 


afAtg _ l / a f\ 

Z S 2 M e (Me-l)\a e j 

or in an approximate form by substituting equations (la) and (2a): 


(39) 


afAtg 



7-1 



Mf 2 


-,-1 

7-i] 



Mf 


(39a) 
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The test time parameter 


of equation (39) is plotted in figure l4. Note 


a f At 2 

l s 2 


that for a given Mf significant gains in At 2 are obtainable by the tunnel 
mode of operation; equation (39 a ) indicates gains roughly proportional to 

Zri 
(A) 2 • 


Also shown in figure l4 is the nonreflected shock tunnel testing time 
together with ticks indicating the values of A required to attain Mf by a 

steady nozzle expansion. The strong-shock approximations require: 



« 7-1 JL 

/ 7 + "l ^f 


(40) 


A, 


NRS 


/ 


2 (7 7 - 1 

7(7 - 1)\ 2 7 + 1 , 


7+1 

, 2 ( 7 - 1 ) 


M f 


7-1 





A NR3 ~ 0-003M.f^ 


(7 = 1-4) 

J 


(4l) 


The length Iq ^ 
required. 


is used for the shock tunnel since there are no 



lengths 


The zero-nozzle-length expansion-tunnel test times are much shorter than 
those for the nonreflected shock tunnel, although the latter, of course, must 
have a much larger area ratio nozzle for the same Mf. The effect of nonzero 

nozzle length will be discussed subsequently. 


Nozzle 

In order to specify the nozzle length ijj, consideration must be given to 

the time lost in starting and stopping the nozzle. An analysis of this problem 
for expansion-tunnel operation will be based on the perfect start approach. 

(See ref 4.) Expressions for the nozzle-time losses will be found for nozzles 
in which the local area ratio varies with both the first and second power of 
distance (linear area ratio and conical nozzles) . 
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The notations used are At flow , At + , and At_ for the times required for 

a transversal of the nozzle hy a fluid particle, a downstream characteristic 
wave (which travels with local velocity u + a) , and an upstream characteristic 
wave (which travels with local velocity u - a) , respectively. The following 
expressions are valid for the times shown in figure 13 : 

For a nozzle where A(g) = A e (l + Cg) 


a f At flow = ^f dg 

° 0 u ( £) 



a(|) M(|) 


(^2) 


where 



(^2a) 


m 


Also 



(kk) 


(^ 5 ) 


a^At + 





l 


W 



(t5a) 
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a^At 

l N 




(46a) 


The nozzle "starting loss" is herein defined as the testing time lost (at 
the nozzle exit) due to the starting process and is equal to the difference 
between the times t x - t w (see fig. 13 ): 


4^ start - ^x ' t w 


At 


start 




t y - t x) + (^y " *vj 


(47) 


^ start “ - ^ 


flow 


^f^ start 


= K- 


N 


(A,Mf) (I 2 + I 3 ) 


(47a) 

(48) 


Similarly, the nozzle "stopping loss" is defined as the test time lost (at the 
nozzle exit) because the downstream characteristic disturbance which signals 
the arrival of the expansion-fan trailing edge at the nozzle entrance precedes 
the fluid through the nozzle: 


4^ stop - ^y' - ^z* 


At 


stop ~ (V _ V) " ( t z' “ t x') 




(49) 


4^ stop ~ 4^tp2. ow - At + 


(49a) 


—f 22 - %(a^) (i 2 - Ij) 


(50) 


The total loss in test time due to the finite nozzle length is then the 
sum of the starting and stopping losses: 


a f At loss 


N 


= 2K- 


( A ?Mf) 


12 


(51) 
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Corresponding equations for a conical nozzle, where 
are as follows: 


A(|) = Ae(l + C'f) 2 , 


where 


a f At flow „ 

= i\s 

% 

a f At + = j 

a-pAt 

— = * 

a f At loss 




A,Mf) (j-L/g " J 5/2) 

(52) 

CVl 

1 

C\J 

rH 

(53) 

! ( J l/2 + J 3/2) 

(54) 

: 2K 2 (A,R f ) J 3 / 2 

(55) 


J n (Me,M f ) 


% [r 2 (M)] 


3-7 

Mr-D 


dM 


(56) 
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M 


K, 




fit 


2 VI - A 2 


[ r a( M f)| 


37-1 

M7-1) 


(57) 


Curves of I n ^l,Mf^ and J n ^l,Mf^ are plotted in figure 15 for y = 1.4 
and 7 = 5/ 3- An arbitrary lower limit for M g of unity has been used in these 
plots; however, the value for any specified limits is simply the difference 

2 


of the values at M = Mf and M = Mg . For the case where 


M T~ 

\ 1 Nj 


» 


7 


approximations to the indefinite integral for I n and J n are 

1 7+1 


I n (M> » (4-i) 


M ?" 1 
7 + 1 
7-1 


- n 


(58) 


- n 
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3-7 

Mr-D 


M^~ n 
7 + 1 

2(7 - 1) “ 1 


( 59 ) 


These approximations are shown as the dashed curves in figure 15- It is evi- 
dent that such approximations would he very good for evaluating the finite 
integral when and M 0 are not small. 

In order to evaluate the fraction of the time At 2 lost in the nozzle 
processes, a nozzle loss parameter is introduced and defined as 

This parameter is plotted for the linear-area nozzle in 
figure 16 together with the approximate equation: 





1 

Zzl 




(60) 



k ^ 


1 


1 - (A) 


3-7 

2 


3-7 M f 


1 - (A) 


(60a) 


Equation (60) retains orders in the reciprocal of A which have been previously 
neglected. This retention is necessary to gain accuracy for this case. 


The nozzle losses are not excessive for 


Z N < l S 2 ‘ 


In fact, for large A, 


the nozzle length can actually exceed the accelerating-chamber length before 
the test time vanished ^that is, Atq oss = At2)- The reason that the 

A = 1 curve does not have At-L QSS = 0 is that by definition the test time is 

bounded by the upstream and downstream characteristics from the nozzle entrance 
(fig. 13)- Consequently, for A = 1 and » 1, the nozzle loss parameter 


approaches two. 


The comparable loss parameter for the nonreflected shock tunnel is also 
plotted in figure 1 6 by using the following approximation: 
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4 yjj + 1 _1_ 

5-77-1 Mf 


(61) 


The values of A of the shock tunnel are indicated by ticks on the curve. 

At a given value of Mp and for the same ratio of %/^gj these curves 
indicate that At^ oss /At 2 is greater for the expansion tunnel than for the 
shock tunnel. However, to obtain high Mp, the value of A becomes extremely 
large in the nonreflected shock tunnel ^for example, at Mp = 50, Ap^g ~ 10 8 j 
so that for a given nozzle wall slope the shock-tunnel nozzle length 

might be much larger than, that of the expansion tunnel. A further considera- 
tion of this problem can be found in the "Discussion" part of this paper. 


The ratio of nozzle loss times for expansion-tunnel nozzles, the radii of 
which vary linearly with distance, to the loss times for those nozzles with a 
linear area variation is for the same Mp, Zjt, and A 


(^l°ss)ii ne ar radius _ K 2^3/2 _ , 

(At los A . V? 

\ iUbb /linear area ^ ^ 


( 62 ) 


where 


F(A) 


2-7 

IT 

5 - 7 1 - (A) 

2(2 - 7 ) -1=2. 

1 - (A) 2 



( 62 a) 


F(A) 


_ -0.3 

M - (A) 

3 ,-.- 0.8 

1 - (A) 


+ (A) 


-0.5 


(7 = 1 . 4 ) ( 62 b) 


For Mp > 20 and A > 10, the approximation values of equation (62b) 
which are plotted in figure 17 are within 3 percent of those of equation ( 62 ), 


V K 1 X 2 ) 


The nozzle losses in the linear-radius nozzle (conical nozzle) are 


thus slightly in excess of the linear-area nozzle and approach a limit of 4/3 
as A approaches <». 
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Another aspect of the nozzle starting process which was considered was the 
nozzle-length requirement to ensure that the accelerating gas flow had heen 
completely established in the nozzle before the test gas reached the nozzle 
exit. This criterion prohibits any interference between the nozzle starting 
processes of the two gases. It was found that this requirement was satisfied 
as long as was no greater than approximately one-half Zg^ for A > 10 

and for values of 7^ equal to both 1.4 and 1.67- 


Intermediate Chamber 

The intermediate chamber length Z Q is determined so that the reflection 

b l 

generated at the intersection of the leading edge of the expansion fan with the 
driver-gas — test-gas interface (between 0) and (2) ) will arrive at the nozzle 
entrance simultaneously with the trailing edge of the expansion fan. (See 
fig. 1.) The governing equation is 


l S± _ 7 + 1 1 


fe)" 

Zg^ 7 - 1 M e - 1 

^1 



3-7 

2(7-1) 


( 63 ) 


is plotted in figure 18. For a given value of A 
decreases rapidly with increasing M (approximately 

as M"5 for large M) . And at a constant (large) value of M, the ratio 
Zgl/Zg^ increases nearly as A®'^. Thus, expansion-tunnel operation requires 

longer intermediate chambers than those for an expansion tube with equal Zg^. 
In the extreme case of large A and low M^, the lengths Z g^ and Z g^ are 
of the same order; however, for high M^. and large A, Zg_^ is only a few 
percent of Zg^. 


and this ratio of Zc 


J 1 


and Zg^, the length Zg^ 


Dump -Tank Length 


An approximation to the dump-tank length Zo may be found by introducing 

R 

the following simplifying assumptions: (a) the nozzle length is zero; (b) the 

shock velocities U a and U Q are constant during shock-wave traversal of 

s 10 S R 

the dump tank; and (c) the total traversal time of these shocks is equal to the 

time interval between the arrival at the nozzle of the shock wave U Q and 

S 11 


the trailing edge of the expansion fan. Thus 
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( 64 ) 




+ At. 


0 - 


l 


s 2 



which can he operated on by substituting equations (3), (28 ) , (32), and ( 39 ) 
plus the relation 


At]_ 



s 2 


to yield 


3 R 


( 7 10 ■ x ) 

( 7 10 + 1 ) 

( 7l1 “ 1 - 1 l 

t 7 n + 1 V 

0 
1— 1 

NA 

- 1 

V 11 + 1 Me - ij 

X + * A 


1/2 


\ 


7 11 " 1 M 


30 


( 65 ) 


The variation of l with is shown in figure 19 (a) for 7^ = ^ 

and 7^ = -k and in figure 19(b) for 7^ = 7^ = — . Larger values of Iq ^ 

are required with increasing A, but the percentage increase is small for 

7 5 

large Mf. Also the values of are smaller when 7^ = j rather than j. 

7 

Not shown are plots for 7 ^q = 7^ = ^ but these values are approximately two- 
thirds those for ^ 7^ = 1 (fig. 19(a)). 


Driver-Chamber Length 

The driver cross-sectional area is assumed to be equal to that of the 
intermediate and expansion chambers. The driver length is determined so that 
the reflection of the driver rarefaction wave off the end plate will pass 
through the point m (fig. l) where the entropy discontinuity between Qj) and (2) 
intersects the leading edge of the main expansion fan. For these conditions 


_ a 4 ^i 
Z Sq U Sl At 0 


( 66 ) 
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(66b) 


From figures 20 which are drawn for values of a 4 ja.^ equal to those of 
figures 12, it is evident that ^dAs^ decreases with increasing A. Conse- 
quently, in terns of matching the driver and intermediate chambers, the expan- 
sion tunnel is less restrictive on driver length. However, since ^Sp/^S 2 

increases with A (fig. l8) , the net result is an increase in ^d/^S 2 A 

34 /l44\ 

as illustrated in figure 21 for the intermediate value of — = 51 )• The 

a f \^9 ) 

driver lengths in the expansion tunnel are thus significantly larger than those 

of the expansion tube; for example, at Mf = 50, 

at A = 100 and approximately 25 at A = 1000. 

Z D « O.Ollg for Mj. = 50. 


N 


ETun 


(M 


is approximately 10 


ET 


However, even at A = 1000, 


DISCUSSION 


Diaphragm Bursting 

Secondary diaphragm .- Reference 1 pointed out the criticalness of the 
secondary-diaphragm burst because of the fact that only an extremely small part 
of the gas initially in the intermediate chamber is actually used for testing. 
If the extent of this gas in state © ^ designated as z, from continuity 


z = — AopAtg 


and by substitution, equation ( 67 ) becomes 


(67) 
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7-1 
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7 + 1 1 
7 - 1 Mf 
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7-1 


M f 
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-1 


7-1 
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7+1 
(A) 2 


(68a) 


Thus, for equal values of Zg 2 and test conditions (f ) , the ratio of the 
length z for the expansion tunnel to that for the expansion tube is 


(z) 


ET 


/ x Z±1 

(Z W» (A) 2 


r l(“f) 


17-1 


7-1 

1 ~ 1 Mf* + (A) 2 


Mj. - 1 


ZzL 

Mf - (A) 2 


(69) 


The relation of equation ( 6 7) is plotted in figure 22 for = 0 and 


N 

shows order-of-magnitude gains in z / l S 2 increasing A. The primary part 

of this gain is attrihutahle to simply the fact that to store the same mass of 
"test slug" in a reduced cross-sectional area requires a longer length if the 
initial charge densities were the same. The other contribution to the net gain 
is a result of the increased test-slug mass due to increased testing time 
more than compensating for the increase with A of the initial charge den- 

7+1 

2 _ 2 

sity p^. For Mg » — , the gain is proportional to (A) (eq. ( 69 )) at 


7 - 1 


a given Mf . 


The secondary-diaphragm bursting problem is alleviated even further when 
diaphragm opening time is considered. If the diaphragm is assumed to he near 
its rupture stress at pressure p^, for a given material the opening time is 

proportional to the tab radius divided by the square root of the primary shock 
pressure ratio P 2 /Pi‘ The fluid mass passing the diaphragm station during 


opening is thus included in a length proportional to 



1 

2 


« a;L d. 


An 


indication of the secondary diaphragm problem is then the ratio z/d with 
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larger values of z/d desired. For the same test-section conditions (f) , test 
section size, and accelerating-chamber length Ig^j the ratios of z/d are 

found from equation ( 69 ) to he 


( z / d )ETun w 

( z A0;ET 


7+2 





(70) 


For A « 100 and A ~ 1000, the ratios of equation (70) are approximately 
2.5 X lO^ and 1.25 X 10^, respectively, and illustrate the large gains in alle- 
viation of the secondary- diaphragm problem for the expansion tunnel. 


The length z test be slightly smaller than 


At. 


z test 


= 'Lr 


test 


At r 


due to the nozzle losses. However, 


Zg ^that is, 

the diaphragm bursting 


affects the downstream part of the test gas, and this gas is also that used in 
the nozzle starting process. Consequently, an imperfect diaphragm burst 
influences first that part of the gas considered expendable in starting the 
nozzle of the expansion tunnel. 


Although the secondary-diaphragm bursting problem can be drastically 
reduced for large A, the problem is not eliminated. If the diaphragm is 
assumed to shatter and the resultant fragments then require acceleration to 
the local free-stream velocity, simple approximations show that the momentum of 
these fragments is not negligible compared with the momentum of the test-gas 
slug (that is, the mass of the test gas slug PpzAg is not orders of magnitude 

greater than the diaphragm mass). An experimental study (ref. 5) using 
poly [ethylene terephthelate/J and cellophane diaphragms for Mg_^ ~ 3 considers 


this opening problem. Since the pressure ratio across the incident and/or 
reflected shock waves increases with Mg_^, the problem for shock speeds perti- 


nent to expansion-tunnel operation is not as severe as at 



3 . If the 


diaphragm petals are restrained so that they fold back against the wall rather 
than proceeding downstream, there is a possibility that the momentum lost from 
the stream may be reduced. Another possibility is the use of electrodynamic 
forces to assist the diaphragm opening. 


Primary diaphragm .- The primary-diaphragm rupture problem should be miti- 
gated by increasing A because of the smaller diameter diaphragm cross section 
If the diaphragm is stressed to rupture at a certain percentage of p^, the 

opening time is proportional to the diameter and will decrease with increasing 
A for given test- section sizes and conditions. 
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Tertiary diaphragm .- A third diaphragm, not present for expansion-tube 
operation, is located at the nozzle entrance of the expansion tunnel. However, 
this diaphragm does not require a perfect rupture since the accelerating gas 
from Zg as well as the initial part of the gas from Zg^ used to start the 

nozzle both pass the diaphragm station prior to the arrival of the gas com- 
prising the test slug. 

Driver Analysis 

The driver pressure ratio Vl^/Vf (fig- 12) and driver length ratio 
Zp/Zg 2 (fig- 21) increase with A. However, these ratios are not a true meas- 
ure of the efficiency of the apparatus since for a given test-section size, the 

_-l/2 

driver cross-sectional area decreases in proportion to A and also for a 

given test time Atg, the length Zgg decreases as A increases. A per- 
tinent parameter which considers all these effects is the energy parameter 

E which is proportional to the energy required in the driver per -unit test- 
section area, test time ( Zjj = 0), and static pressure: 


E = 


E 


ApP^a^Atg 


(71) 


E = 


71). - 1 Pf 


M _ Wd 


T], Z 


hr Sc 


-i 


(j a fAtg ) 

l S 2 


(71a) 


where 


E 


= K 1 )* 1 - |) A e ; 


D 


(72) 


and Tp is the initial temperature before energy addition in the driver. 

— _ >5 

Figure 23 shows the dependency of E on both Mp and A for 7 ^ = j 
Tk Tli 

and values of — = — equal to those of figures 12 and 20. Lower values of 
T i T f 

E result as A increases at a given Mp and T^Tpj and the reductions 
become more pronounced as T l)./Tp increases. More efficient operation appears 
to occur for the higher Tl)./Tp value at the same Mp and A. Thus, when the 
expansion tunnel is compared with the expansion tube on the nondimensional 
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energy basis E, the advantages of the expansion-tunnel mode of operation are 
very significant. This energy reduction is an important consideration for arc 
drivers which use capacitor bank discharges, inductive storage systems, and so 
forth . 


Efficiency 

The ideal efficiency of the expansion tunnel is defined as 


^o = 


Test- section energy 
Driver energy 


_ (| U f 2 )(Pf U f A f At 2) _ y M f 5 


E 


2 E 


(73) 


This efficiency is "ideal" because there are no nozzle starting losses, 
no scoop-off losses, and so forth considered. Figure 2k shows that although_ 
r) Q decreases with Mp for fixed A, it increases very significantly with A 

Th 

for a given Mp. For — = 25 and Mp = 50, t] q equals 0.003 for A = 1, 

T f 

0.025 for A = 100, and 0.053 for A = 1000. 

If the physical dimensions of the nozzle are such that ^ < /a", the nozzle 

entrance section must be designed to capture only that part of the gas in the 
accelerating chamber enclosed within a stream tube of diameter d'. (See 
sketch (a) . ) 
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Consequently, the overall efficiency which is the product of the ideal effi- 
ciency, the nozzle starting efficiency, and the nozzle capture efficiency is 


71 = Vi K V 


7 


= \ 


^test /d 1 
At g \ d / 





m 


Viscous Effects 

Another problem pointed out in reference 1 was the viscous effect in the 
accelerating chamber. One measure of the effect of viscosity is the Reynolds 
number based on diameter. The ratio of NRe,d in the accelerating chamber to 

Nr 6 ,d of the test section (which is also the Reynolds number ratio for the 
accelerating chambers of an expansion tunnel to that of an expansion tube for 
the same Af and Mp) is 


(%e,d) 


accelerating chamber Pe u e^ Pf 


(%e,D) 


test section 


and for p <x T 0 - 1 


(%e,d) 


( N Re,D) 


test section 


He P-pn^D 


accelerating chamber .-2 
_ « (A) 


i- ( 7 -l)o) 


_d_ 

d' 




J 


(75) 


where 


A - 


D_ 

d' 


3 

Equation (75) is plotted in figure 25 for o> = — , which is a typical value for 

4 

air. This figure shows small gains in the Reynolds number ratio with 
increasing A when d' = d. 


Another measure of the effect of viscosity is the ratio of 5/d where 
5 is the boundary- layer thickness which would develop in a length on a 

semi-infinite flat plate with an external flow equal to the flow conditions © 
or (D ; thus , 
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d P e u e d 1 
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(76) 


Thus, an expansion tunnel and an expansion tube which both have the same test- 
section diameter D can be compared by substituting equations (2a), (39 ) > and 
(75) into equation ( 76 ): 

—, 1/2 

•'V r ' 1) (“'*) (At2> 
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and for Mg » 1, Mp » 1, and 1 


/ETun /d 1 
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In contrast to equation (75) the ratio of equation (77a) indicates a 
slight penalty for the expansion tunnel when judged on this basis. (See 

fig. 25 .) Even though the diameter is reduced drastically (<xA ^ ), the unit 
Reynolds number is increased and the length l is reduced sufficiently to 

produce, when d’ = d, only a small change in the 5/d ratio which is propor- 
tional to 

For the limiting case of A » 1, M e » 1, and Mp » 1, equation ( 76 ) 
may be put into the alternate forms: 
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2 Pf 
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f u f 
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( 76 a) 
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fAtg % Mp/d’ \ 
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( 76 b) 


d^ 

For a given condition © and the same values of d, D, and either Zg 

or Atg, significant reductions in 5/d result by reducing d' (increasing 

A). For 7 = 1.4- and co = 3/4-, 5/d is proportional to the 0.7 and 0.9 power 
of d’/d in equations ( 76 a) and ( 76 b), respectively. Any reduction in 5/d 
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resulting from decreasing d’/d is obtained at the expense of reduced nozzle 
capture efficiency and requires an increase in driver pressure and energy. 

This penalty may only be critical at maximum facility performance since excess 
driver capability would be available at "off-design" conditions. One must also 
consider the danger of transition to turbulent flow as d'/d is decreased 

P e u e d 

since the Reynolds number increases. If transition occurred, the asset 

Pe 

of a decreasing laminar ratio of S/d with decreasing d'/d would become a 
turbulent liability. 

Thus in regard to the viscous effects, there do not appear to be either 
large gains or penalties for the expansion tunnel compared with the expansion 
tube for d' = d. However, for d' < d, there is a possible alleviation of the 
viscous effects at the expense of performance capability. 


Low Vacuum Considerations 


Accelerating chamber .- The pressure level Ppp/Pp in the accelerating 

chamber increases with A (fig. 7) and consequently reduces the low pressure 
requirements in this section. This reduction offers significant gains in 
another practical aspect for the expansion tunnel over the expansion tube. The 


low values of Ppq/Pf for A = 1 at high Mf> may restrict the gas in 

state (l^ to helium or hydrogen which have high sound speeds. (See ref. 1.) 
However, for the expansion tunnel, the level of Vii /Vf is raised sufficiently 


so that the same gas can be used both for the test gas in state © and the 
accelerating gas in state (l^ . Thus any diffusion of the accelerating gas 
back into the test gas across the interface between (e) and (2^ would be 
particles of the same molecular species which, through collisions, should soon 
reach the ambient temperature of the test gas. Thus the diffusion contamina- 
tion problem could be curtailed. 


Another advantage of using the same gas in states (T) and ( 0 ) would be 
that the gas in the boundary layer on a model would consist of the correct test 
gas particles even before arrival of the test gas itself. Consequently, the 
boundary layer could probably equilibrate more rapidly to the new test gas flow 
condition than if it initially contained a foreign gas in state (0) . 


Dump tank .- For the boundary conditions of a "perfect start" for the 
accelerating gas, the nozzle charging pressures Ppo/Pf required drop signif- 
icantly with A. (See figs. 8 and 9 and eq. (23).) However, it may be possible 
to raise the pressure above the p-^Q values shown and accept an initial imper- 
fect start because the perfect start conditions produce pressures and velocities 
which generate upstream expansion waves in the dump tank when the accelerating- 
gas — test-gas interface arrives. Thus an increase in p-^Q might be tolerated 

with a reduced- strength dump-tank expansion wave. It is not known whether the 
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pressure could be raised to the point where this expansion wave had zero 
strength ^p g = p^ without slowing down the accelerating-gas nozzle start by an 

unacceptable amount. 

A further consideration of the pressure p is that it does not have to 

meet an exact value but only an upper limit. Thus, in contrast to the expan- 
sion tube where the lowest charging pressure must be closely regulated since 
it determines the strength of the primary expansion wave and hence the test 
conditions, for the expansion tunnel only an upper limit need be set at this 
lowest pressure P]_q - This condition is very desirable from an operating 
viewpoint . 


Effect of Nozzle Configurations 

General effect .- The nozzle design for expansion-tunnel operation is sig- 
nificantly different from that in conventional wind tunnels or shock tunnels 
because the entrance Mach number M e is very large. The boundary condition of 

both large entrance and test-section Mach numbers requires relatively long noz- 
zles. For example, if Mg = 10, the flow along the nozzle center line does not 

even start to expand until the distance downstream of the nozzle entrance is 
equal to 10 accelerating chamber radii. 

If a variation in test- section Mach number is to be obtained simultan- 
eously with no flow inclination or velocity gradient in the test section, con- 
toured nozzles are required for all combinations of Mg and Mp. Each such 

nozzle would have to be corrected for viscous effects which produce both a non- 
uniform entrance velocity profile as well as a nozzle boundary layer. Rather 
than consider such a multiplicity of nozzles, the alternative of conical noz- 
zles is treated briefly in this section with the additional restriction of 
uniform entrance conditions. 

An effective aerodynamic area ratio A is defined as the value of A 
determined from equation (l) with M^ = Mp, where M^ is the nozzle center- 

line Mach number. The geometric area ratio A geom is the ratio of the nozzle 

cross-sectional areas, measured perpendicular to the nozzle axis, at the nozzle 
exit and entrance 



Geometric and fluid mechanic parameters .- At this point three parameters 
are introduced and defined as (see sketch (a) ) 
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(79b) 


These G parameters are a combination of pure geometric factors modified 
to various degrees by aerodynamic factors. The parameter G 0 is principally 

dependent on the geometry of the accelerating chamber Zg^/d and nozzle wall 

angle; whereas the parameters G and G* are more strongly influenced by 
fluid-mechanic aspects. The asterisk superscript denotes critical values when 
the lengths Zg^ and Z^ are such that Atp oss = At 2 (i.e., Atp es p = 0); 

substitution of equation (78) into equation (79) then gives 


G = 




G 



(80) 
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The value of G for the nonreflected shock tunnel with a conical nozzle 


is 
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( 81 ) 


For the nonreflected shock tunnel the reference length is Zg rather than 
Zg^. Values of G* are plotted in figures 26(a) and 26(b). An increase in A 

■X" -)£■ 

requires an increase in G , values of G for the expansion tunnel ranging 
from zero at A = 1 up to approximately 10 at A = 10^. For the nonreflected 
shock tunnel, G* increases as Mpb-5 (eq. (8l)) and varies from 10 at 

= 10 up to approximately 125 at Mp =50. Note that G* is principally a 
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function of A until the value of A becomes large enough to require Mg to 
no longer be much greater than unity. (See fig. 26(b).) 

The testing time for a given configuration is simply 


a f At test 

2q 

b 2 




( 82 ) 


This test time parameter is shown in figure 27 at values of G of 10, 
100, and °o for the expansion tunnel and at these three values of G and 
also at G = 200 for the nonreflected shock tunnel. For A = 103, G = 10 
results in a substantial reduction in test time with increasing Mf, and in 
fact At-^gg-t approaches 0 at Mf « 40. (See fig. 26.) The penalty for 
A = 103 is reduced to only about 10 percent at G ~ 100. (See fig. 27 .) 
When A = 10 2 , the G = 10 case shows approximately 50 percent penalty; and 
for G = 100, the reduction is 5 percent. 


The shock tunnel is more sensitive to low values of G with larger penal- 
ties in test time, the latter vanishing at Mf ~ 9 for G = 10 and Mf « 

for G = 100. 


af Attest — 

The test time parameter depends principally on Mf and A for 

l S 2 

— ax.At +oe+ 

Mf » 1, A » 1, and of order of 1. The relation of — ■ - - to var- 

l S 2 l S 2 

ious other expansion-tube parameters is illustrated by the following set of 
approximate equations. The second term on the right-hand side is the correc- 
tion for the case when Mg 'fcfi 1 . 
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is neglected and Mg » 1, these equations 


If the effect of A / A g e om 


indicate 

0 a. _ — J- 

that Mp 1 will increase when 

^ Zq 

bg 

(a) 

A increases for fixed Zj^Zg^, 

(b) 

^n/^S 2 decreases for fixed A, 

(c) 

G increases for fixed A or Zjj/Zg^, 

(a) 

d'/d decreases, and 

(e) 

9 W increases.. 


These gains in test time would he modified somewhat if consideration were 
given to the relation between A and A g eom j the indicated beneficial effect 

due to increasing B w (eqs. (83e) to (83h)) probably suffering the largest 

reduction. Note for positive test times that the maximum ratio of d'/d 
is specified by equations (83) for given values of Zjj, Zg^, ®w> D, and d. 

It may also be expressed as a function of G 0 and G* (by using eq. (82)) as 



G 0 < G* 


or 


d 

d 


G 0 

= 1 for -2- > 1 
G* 


The parameter G Q is limited by the physical operating conditions of the 
apparatus. Because of attenuation or other various reasons, the factor Zg /d 
will usually be less than 200. Values of tan 9 W will probably be close to 
0.1, although a value of 0.2 might be acceptable for some uses. The third com- 

\fx - 1 

ponent term of G Q , — , is a function of Mg, 
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9 . 


A c 


- 1 


w> 


and the distance 


■'geom 

downstream from the nozzle entrance. This term has a value of zero at the 
entrance and maintains this value until the first characteristic wave from the 
nozzle entrance reaches the center line. The term will then increase in value 
with distance and reach a maximum appreciably in excess of unity, and thereafter 
decreases toward the asymptotic conical nozzle value near unity. Thus a typical 
value for G 0 might be 20. 


The test-section size must also be considered since practical aspects pro- 
hibit the use of very large test-section diameters. Consequently, the range 
of the ratio of D/d will be bounded and, as a result, very large values of A 
can be obtained only at the expense of decreasing d'/d. 
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Growth of the nozzle boundary layer for the rather low Reynolds numbers of 
the nozzle is also a factor for consideration. However, only fairly small 
effective nozzle half-angles are required. Consequently, at first glance the 
problem appears to be principally one of correcting for the effect of d 6 */d£ 
at small £, and then at a larger £ expanding further, if necessary, to 
assure a potential "core" at the test section. This low Reynolds number nozzle 
problem will require additional study before its full effect can be determined. 


Example 

A numerical example, illustrating the discussions in the foregoing para- 

Z q 

/ \ / \ 13 

graphs, is solved with the assumptions: (a) = 150 , (b) — = 10 , 

d a 

(c) G 0 = 20, and (d) A = Ag eom . Equation (82) may then be rewritten as 


V^test = a fAt 2 pS 2 \ / d 
D 2q„ \ d YD 



(84) 


a f At test = 15 a f At 2 L _ 1 Gl_ 
D Z S P \ 2 \Tx 


(84a) 


Equation (84) is in a form which permits evaluation and comparison on a 
test-section-diameter basis which is significant when considering large A. 

a f Attest 


The variation of 


D 


with Mf is shown in figure 28. Since 


d' -- 1 / 2 

— = 10A , no curves are drawn for A less than 100. The gains with 

d 

increasing A of At^. es ^. are larger when compared on a common test-section- 

diameter basis (fig. 28 ) than on a common accelerating -chamber-length basis. 
(See fig. 27-) This statement is not universally valid but applies under the 
previous restrictions (a) to (d) . 

Equation (83e) may be rewritten for Mg » 1 as 
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Examination of equation ( 85 ) reveals the reason that the shock tunnel (operating 

— , afAttest 

at very large A) has a higher value of 


D 


than the expansion tunnel 



r 


for which curves are drawn only for 


A % 10 5 . 


„ , , . „ a f^test „ 

Evaluation of from 

D 


equation (84), with the corresponding A(M) of the shock tunnel, will not pro- 
duce an answer in agreement with the shock- tunnel curve for two reasons. The 
first is that the assumption Me » 1 is violated when M 0 approaches Mg,> 

and this violation introduces an error which results in the first term of equa- 
tion (84) being low by a factor of 5-8. Even if this correction factor were 
taken into account, the expansion tunnel with Mg = M 2 would not have the same 

afAto — 

value of as the shock tunnel at the same Mp and A. Instead 

D ^ 


a^Atg 
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would be greater than 
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by the factor of 6.74 which is 
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This variance should not be interpreted to mean Atgrj-^ > 


because 
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D / ETun 


was assumed to be equal to 
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whereas the overall 


NRS 


length Iq ^ + l of the expansion tunnel is 7-74Zg^. Consequently, if a 


com- 


/ETun 


would be equal to 


. 87(At 2 ) NR g- This answer is different from unity only because by the 


pari son were to be made on a total length basis, (Atp^ 

. 0 . 

i.ju 

definition of (At 2 )E>p un ^ gas i- n condition (2^ is not included; yet for 
Me = Mg,, state (g) = state (l^ and state (g) = state (e) = state (2^ . Con- 
sequently, if this gas in state (3) were also used, the times for the two 
apparatus would then be equal. 


The overall efficiency for this particular example may be expressed as 
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1 G* \l00 
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(86) 


At Mp = 50, the values of T) are 0.019> 0.0070, and 0.0044 for A = 100, 

500, and 1000, respectively. Thus the nozzle capture losses more than offset 
the gains in T] 0 (fig. 24) with increasing A, and the most efficient oper- 
ating point for this example is the value of A such that d' = d (A = 100) . 
Note that the expansion tube (A = l) has an ideal efficiency of only 0.003 for 
M f =50. 


Since t] was not computed for the nonreflected shock tunnel, the exact 
value of T] cannot be found. However, since the nozzle capture efficiency at 
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it is obvious that the overall effi- 


Mf = 50 


is approximately 


100 

A 


10-^ 




ciency is well below the expans ion -tunnel values. 


Design Details 


The design of the conical nozzles is an important problem area where a 
study that uses the method of characteristics has been underway for an extended 
period. The general results of this study are not included here but one par- 
ticular case is illustrated in figures 29 and 30. The entrance Mach number is 
10 and the primary nozzle angle 9 W is 5-3° • In an effort to reduce the noz- 
zle lengths Zjj, the use of compound or multiple nozzles has also been con- 
sidered. This latter scheme, similar to that used in shock tunnels (refs. 2 
and 6), is shown in sketch (b). The curves of figure 29 show the center-line 

Mach number distribution as a function of Zjj/d 
for both the simple nozzle and the compound 
nozzle with the imbedded nozzle lips located 
£■ 

at — = 12. This imbedded nozzle has an 
d 




for a compound nozzle and at 


entrance radius equal to 0.37<I' and a wall 
angle of 5-6° The compound nozzle arrange- 
ment is successful in shortening the nozzle for 


a given ^for example, M^ = 40 at = 40 

1 \ 

— = 75g- for a simple nozzle J. However, this 


reduction in Zjj is obtained at the expense of test-section gradients as 
shown by the radial distribution of flow angle and Mach number in figure 30. 

The velocity variation is negligible since the Mach number 
variation is almost entirely due to changes in the speed 
of sound. Values of 9 for source flow with virtual ori- 
gins at the nozzle entrances indicate that the flow inclina- 
tion is slightly worse than that of source flow. For radi- 
ation experiments where the important contributions arise 
at the model surfaces most normal to the free stream, the 
variation in 9 would probably not be important. 

A practical nozzle configuration might also include a 
boundary-layer scoop at the entrance section to obtain a 
more uniform entering velocity profile. (See sketch (c).) The trade-off in 
nozzle performance against this increased complexity has not been explored. 




Critical Lengths 

Accelerating-chamber critical length . - Since operation over a range of Mp 
and A is desirable and the values of d and D are generally fixed, the 




nozzle capture diameter d* will probably be the variable parameter of an 
actual expansion tunnel. This variation is easily accomplished; for example, 
a reduction in d' results from a slight forward extension to the nozzle sec- 
tion. Thus the "constant" component dimensions are and D. Expressions 

relating the critical acceleration-chamber length Z g for which the test time 

vanishes to these components may be obtained from equations ( 78 ) , ( 79 ) > 
and ( 80 ) . Thus , 
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From figure 26 at IVh = 50, ~~ 0.J2 and 0.77 for A = 10-^ and 

A A 

— p ^S 2 ✓a 1 

A = 10 , respectively. Consequently, for these cases ~ 0.64 and 0-94; and 

7 * l N 

t* a 

3-2 and 7.7 (for tan 9 W = 0.1). The critical accelerating -chamber length 
is not excessive; therefore, practical lengths for Zg^ do not result in large 

percentage of losses in test time. If in the cases above Zg ~ 150d ~ 15D, 

7* 

1 s 2 3 _ ^ 

the test-time percentage loss is ~ — ~ 20 percent for A = 10-'. 

*S 2 15 

Dump-tank length ■ - From equation ( 65 ) and figures 19(a) and 19 (b), it is 

increases with increasing K. Such an increase 


evident that the ratio 


\/ ls 2 


is not a true penalty for expansion-tunnel operation because of the increased 
testing time available. For example, by combining equations ( 65 ) , (39 ) > (3a), 
and (la), the following relation may be found: 
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Consequently, on a unit test time basis, Zg^ actually decreases with A. Of 

course, to utilize fully the expansion-tunnel capabilities the ratios for 

Zg /Zg of equation ( 65 ) should be maintained. 

R / 2 
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RESUME 


The perfect-gas analysis of this report has covered many facets of 
expan si on -tunnel operation. The more important advantages and disadvantages 
of the expansion tunnel compared with the expansion tube are as follows: 


Advantages: 

(l) The zero-nozzle-length test time ^a^Atg/^gg^ P er un; ^ length is 

increased. The time lost in the starting and stopping processes of the nozzle 
is, in general, only a small fraction of this time Atp for ratios of nozzle 

length to accelerating-chamber length Zgg l ess than l/2 and for large 

area ratios A. 


(2) The usable test slug length prior to secondary-diaphragm rupture 

7+1 


increases rapidly with A and is roughly proportional to A 
section Mach number Mp. 


at high test- 


(3) The secondary-diaphragm bursting problem is greatly reduced. 

(4) The nondimen sional energy parameter Ey/ApPpapAt 2 decreases with 

increasing A. Thus the expansion tunnel is especially suited to arc-heated 
drivers. The ideal efficiency also increases with A. 

(5) The primary diaphragm is of smaller diameter for the same test-section 
area Ap and consequently this bursting problem is also reduced. 

(6) The low pressure in the accelerating chamber p^ is significantly 

higher. As a consequence, either the pumping capacity could be reduced and/or 
the same gas as the test gas might be used to reduce interface mixing effects. 

(7) The dump-tank length parameter Zg R /apAt 2 is reduced. 

(8) Viscous effects may be reduced for nozzle entrance diameter less than 
accelerating-chamber diameter (d 1 < d) . 


Disadvantages : 

(1) Fairly long nozzles are required. The (theoretical) expansion tube 
flexibility of variable Mach number Mp without nozzle changes is lost. How- 
ever, conical nozzles might be employed to regain this flexibility with small 
flow gradients. The effect of nozzle boundary-layer growth must also be 
considered . 

(2) The ratio of driver pressure to test-section pressure p^/pp is 
increased . 


I 

(3) The ratio of the minimum charging pressure in the dump tank to the 
test-section pressure Ppo/^f ^ ecrease< ^ an ^ thus requires increased pumping 

capacity. 

(4) A third diaphragm is added to apparatus. 

(5) The length parameters Z^Z Sq’ an< ^ t ncrease 

with A as well as Zg^/a^^zM^ and ^p/^Atg where Z D is the driver length, 
Zg , the intermediate -chamber length, Zg^, the dump-tank length, and a f is 
speed of sound in test section. However, for large M^., the length Zg^ is 
still by far the predominant length as long as A > 10^. 

In the opinion of the authors, the advantages of the expansion tunnel 
appear, particularly in the light of practical operating problems, to outweigh 
the disadvantage s . 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 2, 1965- 
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TABLE I.- INDEX TO FIGURES 


Schematic and wave diagram for expansion tunnel . 

Variation of conditions at nozzle' entrance with test-section Mach number 

( a ) Speed- of- s ound rat io • 

(b) Flow- velocity ratio. 

(c) Static-pressure ratio. 

(d) Nozzle -entrance Mach number. 

Variation of conditions in region (2) with test- section Mach number 

(a) Speed-of- sound ratio. 

(b) Flow- velocity ratio. 

(c) Static-pressure ratio. 

(d) Reciprocal total- enthalpy ratio. 

Variation of parameter |3 with Mach number and area ratio 

(a) Mach number. 

(b) Area ratio. 

Variation of primary shock velocity parameter with test- section 

Mach number 

Variation of initial intermediate-chamber density ratio with test- section Mach number 

a l 49 5 

Initial accelerating-chamber pressure ratio parameter for = — 7, n = — ... 

a ll lW 11 3 

Dump-tank charging pressure parameter for perfect nozzle start 

( a ) 7 10 = 7u = 1-4. 

(b) 7 10 = |j 7u = 1-4. 

(c) 7 10 = 7 X1 = j- 

Dump-tank absolute charging pressure for perfect nozzle start 

Unsteady waves generated by entropy discontinuity flowing through nozzle 

(a) Schematic of wave system. 

(b) Wave diagram. 

(c) Pressure-velocity diagram. 

Pressure ratio across waves generated by passage of entropy discontinuity 

between (e) and (21) through nozzle 

(a) 7ll = 1.4. 

(b) r u = |. 
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TABLE I.- INDEX TO FIGURES - Concluded 


Effect of nozzle length on loss of test time for linear area nozzle ..... 

Ratio of nozzle time losses in linear-radius nozzle compared with those in linear-area 

nozzle 

Ratio of intermediate -chamber length to accelerating-chamber length 

Ratio of dump -tank length to accelerating-chamber length ■ 

7 io = f ; 7 u = 5' 

(b) 7 10 = 7 X1 = |- 

5 

Ratio of driver length to intermediate -chamber length for helium driver. 7^ = — 

(a) ^ = 10. 

Tf 

(b) | = 25 . 

( = ) ^ = 50 . 

i f 

5 -^k 

Ratio of driver length to accelerating-chamber length for helium driver, y = — • — — = 25 . . 

4 3 T f 

Ratio of initial length of test gas slug to accelerating-chamber length. Zjj = 0 

Variation of driver energy parameter with test-section Mach number for helium-driven 

expansion tunnels . 

(a) ^ = io. 

(b) = 25. 

T f 

(c) L = 50. 

Tl 

Ideal efficiency of helium-driven expansion tunnels. — — = 25 

•^f 

Ratios of Reynolds numbers and nondiraensional boundary-layer thicknesses in accelerating and 

test chambers for case d' = d 

Variation of values of parameter G* for conical nozzles ’. 

(a) With Mach number. 

(b) With area ratio. 

a„At, . 

Effect of G on test time parameter — — x,es1: 

l S 

afAt+. p qi. 

Variation of test time parameter — — with test-section Mach 

D 

l S 

number for numerical example. Conical nozzle; ^ = 150; and G 0 = 20 

Mach number variation along center line of simple conical nozzle and compound conical 

nozzle with distance from nozzle entrance. Mg = 10; 0 W = 5.3° 

Radial variation of flow Mach number and inclination 
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Time 
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Driver | chamber 1 

Accelerating chamber 
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I Nozzle 


Shock wave 

Particle path 

Characteristics 
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wave / 


/ / 
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Figure I.- Schematic and wave diagram for expansion tunnel. 






Equation (2) 


Equation (2a) 



10 20 30 40 5( 

Test-section Mach number, 

(a) Speed-of-sound ratio. 


Variation of conditions at nozzle entrance with test-section Mach number. 


Ratio of entrance flow velocity to exit flow velocity 












Test -section Mach number, M-j 


(b) Flow-velocity ratio. 
Figure 3.- Continued. 




Static-pressure ratio, zA 



(c) Static-pressure ratio. 
Figure 3.- Continued. 



Reciprocal total-enthalpy ratio, -==- 



Test -section Mach number, M| 


(d) Reciprocal total-enthalpy ratio. 
Figure 3.- Concluded. 
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Density ratio 



0 10 20 30 40 


Test-section Mach number, Mf 

Figure 6.- Variation of initial intermediate-chamber density ratio with test-section Mach number. 












Absolute pressure parameter, -g— 



Test-section Mach number, M-f 
Figure 9. - Dump-tank absolute charging pressure for perfect nozzle start. 
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Driver pressure ratio, — - 




Figure 12.- Continued. 








7 = 5/3 


7= 7/5 



Test-section Mach number, Mf 


(a) Ij. 

Figure 15.- Integrals I n and J n for y = | and | with lower limit M e = 1. 













Dump-tank length ratio, 
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Figure 19. - Ratio of dump-tank length to accelerating-chamber length. 


84 




Driver length ratio 










Driver length ratio 





Driver energy parameter, 






Driver energy parameter, 




Figures.- Concluded. 


Ideal efficiency, 
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Figure 25. - Ratios of Reynolds numbers and nondimensional boundary-layer thicknesses in 



(a) With Mach number. 

Figure 26.- Variation of values of parameter G* for conical nozzles. 
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Test time parameters, 





Test time parameter 



Test-section Mach number, Mf 


a,Ul. . 

Figure 28.- Variation of test time parameter — with test-section Mach number for 

d 1 \ 

numerical example. Conical nozzle; "q = jq; -jj - = 150; G 0 = 20. 
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Figure 29.- Mach number variation along center line of simple coni 







Figure 30.- Radial variation of flow Mach number and inclination. 
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